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Abstract. Let £ be a fixed vertical lattice line of the unit triangular lattice in the plane, and 
let Ti be the half plane to the left of I. We consider lozenge tilings of TL that have a triangular 
gap of side-length two and in which i is a free boundary — i.e., tiles are allowed to protrude 
out half-way across £. We prove that the correlation function of this gap near the free boundary 
has asymptotics j^p, T — * oo, where r is the distance from the gap to the free boundary. This 
• parallels the electrostatic phenomenon by which the field of an electric charge near a conductor 

£SJ , can be obtained by the method of images. 

o 

(N 

On . 

1. Introduction 

The study of the interaction of gaps in dimer coverings was introduced in the literature 
by Fisher and Stephenson [13]. This pioneering work contains three different types of gap 
interaction in dimer systems on the square lattice: (i) interaction of two dimer-gaps (equiv- 
alently, interaction of two fixed dimers required to be contained in the dimer coverings); (ii) 
interaction of two non-dimer-gaps (specifically, two monomers), and (Hi) the interaction of 
a dimer-gap with a constrained boundary (edge or corner). 

The first of these types of interactions was later generalized by Kenyon [18] to an arbitrary 
number of dimer-gaps on the square and hexagonal lattices, and recently by Kenyon, Ok- 
ounkov and Sheffield [20] to general planar bipartite lattices. Interactions of the second type 
were studied by the first author of the present paper in [4] [5] [6] [8] [7] , where close analogies 
to two dimensional electrostatics were established. 
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Two instances of interaction of non-dimer-gaps with constrained boundaries can be found 
in [19, Section 7.5] (interaction of a monomer with a constrained straight line boundary on 
the square lattice), and respectively [5, Theorem 2.2] (interaction of a family of triangular 
gaps with a constrained straight line boundary on the hexagonal lattice). 

In this paper we determine the interaction of a triangular gap with a free straight line 
boundary (i.e., dimers are allowed to protrude out across it) on the hexagonal lattice. This 
type of interaction seems not to have been treated before in the literature. (We are aware 
of one other paper, namely [10], addressing the asymptotic behavior of lozenge tilings under 
the presence of a free boundary, but the regions considered there contain no gaps.) We find 
that the gap is attracted to the free boundary in precise analogy to the (two dimensional) 
electrostatic phenomenon in which an electric charge is attracted by a straight line conductor 
when placed near it. 

This develops further the analogy between dimer systems with gaps and electrostatics 
that the first author has described in [5] [6] [8] [7] . More generally, our result shows that in 
any physical system that can be modeled by dimer coverings, a gap will tend to be attracted 
to an interface corresponding to a free boundary. This effect, purely entropic in origin, 
is reminiscent of the Cheerios effect by which an air bubble at the surface of a liquid in 
a container is attracted to the walls [32] (note that the Cheerios effect is not entropic in 
origin). 

2. Set-up and results 

There seem to be no methods in the literature for finding the interaction of a gap "in a 
sea of dimers" with a free boundary. However, as V. I. Arnold said, "mathematics is a part 
of physics where experiments are cheap." We now design such an experiment in order to 
determine the interaction of a gap in a dimer system on the hexagonal lattice with a free 
boundary. 

Consider the tiling of the plane by unit triangles, drawn so that one family of lattice lines 
is vertical. Clearly, the hexagonal lattice can be viewed as the graph whose vertices are the 
unit triangles, and whose edges connect precisely those unit triangles that share an edge. 
Dimers on the hexagonal lattice are then (unit) lozenges (i.e., unit rhombi) consisting of 
pairs of adjacent unit triangles. 

The free boundary we choose is a lattice line I — say vertical — on the triangular lattice, 
to the left of which the plane is covered completely and without overlapping by lozenges, 
except for a gap <2 in the shape of a triangle of side- length 2, pointing to the left; the 
lozenges are allowed to protrude halfway across the free boundary, to its right (Figure 4 
pictures a portion of such a tiling; the dotted lines should be ignored at this point) . 

We define the correlation function (or simply correlation) of the hole <2 with the free 
boundary I as follows. Choose a rectangular system of coordinates in which I is the y- 
axis, the origin is some lattice point on £, and the unit is the lattice spacing. Let <2(k) 
be the placement of <2 so that the center C of its right side has coordinates (— k\/3 : 0) 
(i.e., C and the origin are the endpoints of a string of k contiguous horizontal lozenges; 
Figure 4 illustrates <2(2), the origin being denoted by O there). Let H n x be the lattice 
hexagon of side-lengths 2n, 2n, 2x, 2ro, 2n, 2x (in counter-clockwise order, starting with the 
southwestern side) centered at the origin (thus H n ^ x is vertically symmetric about £, and 
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its horizontal symmetry axis cuts <2(k) into two equal parts; for example, the boundary of 
the region in Figure 3 is #4,4). Let F njX be the region obtained from the left half of H UjX 
by regarding its boundary along I as free (i.e., lozenges in a tiling of F njX are allowed to 
protrude outward across £). Figure 1 shows the region i 7 ^ together with such a lozenge 
tiling; the origin is labelled by O. 




Figure 1 



Following in the spirit of [13] and [5], for any fixed integer k > we define the correlation 
of <2(k) with the free boundary £, denoted u/(k), by 

„,(*) := lin. M(F ^f\ (2.1) 

where M(R) stands for the number of lozenge tilings of the region R (if R has portions of 
the boundary that are free — as in our case — then it is understood that what we count is 
tilings in which lozenges are allowed to protrude out across the free portions). A tiling of 
^4,4 \ <2(2) of this type is illustrated in Figure 4. 

We note that, by [9], lozenges have maximum entropy statistics (in the scaling limit) 
at the center of a regular hexagon. According to this, (2.1) is a natural definition for the 
correlation function. An analogous definition was used in [5]. 

In Lemma 13 (with £ = 1) we obtain an exact expression for Wf(k) in terms of an integral. 
What affords this is an exact formula for M(F nyX \ <2(k)), which we present in Theorem 4. 
We then deduce the asymptotics of uJf(k) as k — > 00 using Laplace's method (see Theorem 15 
with £ = 1). The result is the following. 
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Theorem 1. As k — > oo, we have 



where d zs i/ie Euclidean distance. 

Remark 1. In fact, our results allow us to determine the correlation of <2(k) with the free 
boundary I in a more general situation, namely when the sides In and 2x of F HjX grow to 
infinity so that x/n approaches a positive real number £ not necessarily equal to 1. This leads 
to the correlation u;/(fc;£) defined in (7.11). It turns out (see Theorem 15) that, for £ ^ 1, 
the behavior of this correlation is exponential in fc. This is in contrast to the behavior of the 
correlation of lozenges on the ("infinitely large") torus, in which case Kenyon, Okounkov 
and Sheffield have shown in [20, Sec. 4.4] that the correlation decays polynomially. More 
precisely, it is shown in [20] that, for dimer models on doubly periodic bipartite planar 
graphs, there can only occur three different "phases" characterized by the behavior of edge- 
edge correlations: liquid, gaseous, and frozen. Our situation is readily seen to be in the 
liquid phase. The liquid phase is shown in [20] to be characterized by a polynomial decay 
of edge-edge correlations, and the arguments there imply that the correlation of holes of 
side 2 also have polynomial asymptotic behavior. It is this fact that is in contrast with the 
exponential interaction of Theorem 15. In fact, if one would define analogous generalized 
correlations for the lozenge tiling models in [4] [5] [6] , one would observe the same phenomenon 
of exponential behavior of correlation if the correlation is not "measured" inside a patch of 
the lattice region in which all three types of lozenges occur with equal probability. This 
hints at the limitation of the transfer of properties of dimer models on the torus to dimer 
models of bounded regions in the plane, which is one of the main driving forces in [20] . We 
plan to address this phenomenon in more detail in forthcoming publications. 

In [7] the first author described how a distribution of fixed holes on the triangular lattice 
defines in a natural way two vector fields. The F-field is a discrete vector field defined at 
the center of each left-pointing unit triangle e, and equal to the expected orientation of the 
lozenge covering e (under the uniform measure on the set of tilings). To define the T-field, 
one introduces an extra "test-hole" t and measures the relative change in the correlation 
function under small displacements of it, as the other holes are kept fixed. One can prove 
(details will appear elsewhere) that in the scaling limit of the lattice spacing approaching 
zero, this relative change is given by the scalar product of the displacement vector with a 
certain vector T(z), where z is the point to which the test hole t shrinks when the lattice 
spacing approaches zero. This defines the second field. 

When these fields are generated by lozenge tilings that cover the entire plane with the 
exception of a finite collection of fixed-size holes (the case treated in [7] and [8]), both the 
T-field and the scaling limit of the F-field turn out to be equal, up to a constant multiple, to 
the electrostatic field of the two dimensional physical system obtained by viewing the holes 
as electrical charges. 

But what if we do not tile the entire plane, but only the half-plane to the left of the free 
boundary £, and we have no holes? 

The above definitions for the F-field and T-field would still work, provided (i) the scaling 
limit of the discrete field defining F exists, and (ii) the scaling limit of the relative changes 
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in the correlation function under small displacements of a test hole exists and is given by 
taking scalar products of the displacement vector with the vectors of a certain field. 

Our exact determination of Uf (k) (see Lemma 13) allows us to verify (ii) for displacements 
along the horizontal direction. <2(k) plays now the role of a test charge. We obtain the 
following result. 

Theorem 2. We have 

Remark 2. By symmetry, displacements of <2(k) parallel to I leave ujf unchanged, so the 
relative change in uif corresponding to such displacements is zero. Thus, provided the field 
T exists, it follows from Theorem 2 that its value at z is 

T « - -2d(br (2 ' 4) 

where i is the unit vector in the positive direction of the x-axis (the 2 at the denominator 
comes from the fact that T arises from the expression on the left hand side of (2.3) divided 
by the product of the displacement, a/3 in this case, and the "charge" of the hole <2{k), 
which is 2; see [7] for details). Note that by [7] we would obtain (up to a multiplicative 
constant of 2) the same field T at z if we look at tilings of the entire plane, with the mirror 
image of our test-hole <2(k) being a fixed hole. This is analogous to the phenomenon in 
electrostatics by which the field created by an electric charge placed near a conductor can 
be obtained by the method of images (see e.g. [11, Chapter 6]). 

The F-field could be determined by an "experiment" analogous to the one we described 
at the beginning of this section. What one needs now is the number of lozenge tilings of 
F UjX \ L(k) : where L(k) is the horizontal lozenge contained in <2(k). This turns out to be 
given by a formula very similar to (3.2), namely by 

w/e, \ rfuw f n + k\ 2 (n - k - 1)1 (x + n - k) 2k +i yi ( 2x + 2s) 4n - 4s+ i 



E — 

J—^ tn+k- 



£o ( n 2 + k k +if ^-(n-k-i - !) !2 (n + k-i + l) n _ k (n + k - i + l) t (2n - i + \) t 
■ ( (x)i (x + i + 1) 

n—k — i — l (x -\- fl ~\- k ~\- V) n —k 

- (x) n - k (x + n + k + l) n _ fc _j_i (x + 2n - i + (2.5) 

In the same way as we derive Theorem 1 from Theorem 4 via Lemmas 12-14, Laplace's 
method can be used to deduce from (2.5) the following result. 



Corollary 3. Let e(k) be the leftmost left-pointing unit triangle of<2(k). Then 



F(e(fc)) ~ o ( ] ) , k^oo. (2.6) 
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Remark 3. Equations (2.4) and (2.6) imply that, in sharp contrast to the case of lozenge 
tilings of the plane with a finite number of fixed size holes, for the half-plane with free 
boundary the fields T and F have radically different behavior: while the former behaves as 
the electrostatic field near a conductor, the latter is zero in the scaling limit. 

Our approach to proving Theorems 1 and 2 consists of solving first the counting prob- 
lem exactly, see Theorem 4. This result generalizes Andrews' theorem [1] (which proved 
MacMahon's conjecture on symmetric plane partitions) in the case q=l. Its proof is given 
in Sections 4 and 5, with some auxiliary results proved separately in Section 6. It is based 
on the "exhaustion/identification of factors" method described in [23, Sec. 2.4]. Finally, in 
Section 7, we perform the asymptotic calculations needed to derive Theorems 1 and 2 from 
the exact counting results. 

3. An exact tiling enumeration formula 

Tilings of the region F njX are clearly equivalent to tilings of the hexagon H n x that are 
invariant under reflection across its symmetry axis I. Counting such tilings was a problem 
considered (in the equivalent form of symmetric plane partitions) by MacMahon in the early 
twentieth century (see [26, p. 270]). MacMahon conjectured that the number of vertically 
symmetric lozenge tilings of a hexagon with side lengths 2n, 2n, 2x, 2n, 2n, 2x is equal to 

( x + 2)271 TT ( 2x + 2g )4n-4s+i -n 
11 (o. a \. ~ ' ^- iJ 



4n— 4s+l 



where (a) m is the Pochhammer symbol, defined by (a) m := a(a + 1) • • • (a + m — 1) for 
m > 1, and (a)o := 1. This was first proved by Andrews [1]. Other proofs, and refinements, 
were later found by e.g. Gordon [15], Macdonald [25, pp. 83-85], Proctor [29, Prop. 7.3], 
Fischer [12], and the second author of the present paper [21, Theorem 13]. 

Our "experiment" - - counting M(i ? n x \ <2(^)) ~ ~ is by the same token equivalent to 
counting vertically symmetric lozenge tilings of H UjX with two missing triangles (compare 
Figures 3 and 4). This is in fact a generalization of MacMahon's symmetric plane partitions 
problem (see Remark 4). 

The key result that allows deducing Theorems 1 and 2 is the following. 

Theorem 4. For all positive integers n,x and nonnegative integers k < n — 1, we have 



E 



, x (k)) _{4k + l\ (n + k)l fr (2x + 2s) 4n ^ s +1 
M(F n>x -{ 2k ) {x + n-k) 2k+1 l\ (2,) 4n _ 4s+1 

n — k — l /!-, 

\2~)i 



i\ (n - k - i - l)! 2 (n + k - i + l) n _ fc (n + k - i + 1), (2n - i + \), 
■ { i {x)i (x + i + l) n _ fc _j_i (x + n + k + l) n -k 

- (x) n -k {x + n + k + l) n _ fc -i-i (x + 2n - i + 1)^. (3.2) 



Remark 4- Replacing x by x — 1, n by n + 1, and k by n, one can see that the above 
formula specializes to MacMahon's formula (3.1). More precisely, because of forced lozenges 



THE INTERACTION OF A GAP WITH A FREE BOUNDARY IN A DIMER SYSTEM 7 




Forced lozenges when the hole touches the left border 
Figure 2 

(see Figure 2), the enumeration problem in the statement of Theorem 4 reduces to the 
problem of enumerating vertically symmetric lozenge tilings of a hexagon with side lengths 
2n, 2n, 2x, 2n, 2n, 2x. 

The proof of Theorem 4 is given in the next two sections. In Section 4, we show that 
M(F nj:B \<2(fc)) can be expressed in terms of a certain Pfaffian. This Pfaffian is then evaluated 
in Section 5. 

4. Lozenge tilings and nonintersecting lattice paths 

The purpose of this section is to find a manageable expression for M(F n ^ x \ <2(k)) (see 
Lemma 6 at the end of this section). In this context, we will find it more convenient to think 
of the tilings of F n>x \ <2(k) directly as tilings of a half hexagon with an open boundary (cf. 
Figure 4) as opposed to symmetric tilings of a hexagon with two holes (cf. Figure 3). There 
is a well known bijection between lozenge tilings of lattice regions and families of "paths of 
lozenges" (see Figure 4), which in turn are equivalent to families of non-intersecting lattice 
paths (see Figure 5). Its application to our situation is illustrated in Figures 4 and 5. The 
origin of the system of coordinates indicated in Figure 5 corresponds to the point O' in 
Figure 4 (note that the bottommost path of lozenges in Figure 4 is empty for the illustrated 
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A symmetric lozenge tiling of the hexagon H n)X with two holes. 

Figure 3 

tiling; the corresponding lattice path in Figure 5 has no steps). 

By this bijection, lozenge tilings of F n>x \< 2 (k) are seen to be equinumerous with families 
(Pi, P 2 , . . . , P2n) of non-intersecting lattice paths consisting of unit horizontal and vertical 
steps in the positive direction, where Pi runs from Ai = (— i, i) to some point from the set 
I U {S 1 ,S 2 }, i = 1, 2, . . . , 2n, with 

I = {(-l,s) : s = l,2,...,2x + 2n}, 

51 = (-2k -l,x + n + k), 

5 2 = (-2k-2,x + n + k + l), (4.1) 

and the additional condition that Si and S 2 must be ending points of some paths. 

At this point, we need a slight extension of Theorem 3.2 in [31] (which is, in fact, derivable 
from the minor summation formula of Ishikawa and Wakayama [17, Theorem 2]). The reader 
should recall that the Pfaffian of a skew-symmetric 2n x 2n matrix A can be defined by (see 
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A lozenge tiling of the region F n>x \ <2(k); the right boundary is free. The 
dotted lines mark paths of lozenges. They determine the tiling uniquely. 

Figure 4 



e.g. [31, p. 102]) 

PM:= J2 s S n7r II A ^ ( 4 ' 2 ) 

n€M[l,...,2n] i<j 

i,j matched in tt 

where M[l, 2, . . . , 2n] denotes the set of all perfect matchings (1-factors) of (the complete 
graph on) {1, 2, . . . , 2n}, and where sgn7r = (— l) "^ 71 "), with cr(n) denoting the number of 
"crossings" of n, that is, the number of quadruples i < j < k < I such that, under it, i is 
paired with k, and j is paired with /. It is a well-known fact (see e.g. [31, Prop. 2.2]) that 

(Pf A) 2 = detA. (4.3) 

Theorem 5. Let {Ai, A 2 , . . . , A p , Si, S 2 , ■ ■ ■ , S q } and I = {I±, I 2 , ■ ■ ■} be finite sets of lattice 
points in the integer lattice Z 2 . Then 

^) = (-i) (I) E( s g n7r )-^ nonint ( A ^ su/ )' ( 4 - 4 ) 
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The paths of lozenges of Figure 4 drawn as non-intersecting lattice paths on I?. 

Figure 5 

where = (A^^i), A^), . . . , A n ^), and 7 ?nonmt (A 7r — > S U I) is the number of families 
(Pi, Pi-, ■ ■ ■ , P p ) of non-intersecting lattice paths consisting of unit horizontal and vertical 
steps in the positive direction, with Pk running from A n ^ to Sk, for k = 1,2, ... ,q, and to 
Ij k , for k = q + 1, q + 2, . . . ,p, the indices being required to satisfy j q +i < j q +2 < • • • < j P - 
The matrix Q = (Qi,j)i<i,j<p is defined by 

Qi,i= (V(A i ^I s )-V(A j ^I t )-V(A j ^I 8 )-V(A i ^I t )), (4.5) 

l<s<t 

where V(A — > E) denotes the number of lattice paths from A to E, and the matrix H = 
(Hi,j)i<i<p, i<j<q by 

H iJ = V(A i ^S j ). 

In the special case when the starting and ending points satisfy a certain compatibility 
condition (called ^-compatibility in [31]), the only permutation tt which contributes to the 
right-hand side of (4.4) is the identity permutation, and (4.4) reduces to [31, Theorem 3.2]. 
In our context, the compatibility condition is not satisfied. However, the same arguments 
that prove [31, Theorem 3.2] can be used to obtain (4.4). (Alternatively, one could use the 
minor summation formula of Ishikawa and Wakayama [17, Theorem 2]. In it, choose m = p, 
r = q, and the skew-symmetric matrix B to be Bij = 1 for i < j — which makes all principal 
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Pfaffian minors of B equal 1 — to expand the Pfaffian on the left-hand side of (4.4) into a 
sum of minors of a certain matrix. Each minor can then be seen to count certain families of 
nonintersecting lattice paths by the general form of the Lindstrom-Gessel-Viennot theorem 
[24, Lemma 1], [14, Theorem 1], and, altogether, these are the families that are described in 
the statement of Theorem 5.) 

We now apply Theorem 5 to our situation, that is, p = 2n, q = 2, Ai = for 
i = 1,2,..., 2n, and Si, S2, and I are given by (4.1). It is not difficult to convince oneself 
that, for this choice of starting and ending points, all families of nonintersecting lattice 
paths counted on the right-hand side of (4.4) give rise to even permutations it. Hence, the 
right-hand side of (4.4) counts indeed the families of nonintersecting lattice paths that we 
need to count. By Theorem 5, their number is equal to the negative value of the Pfaffian of 

M n (aO:=(_£ *y (4.6) 

where Q is a (2n) x (2n) skew-symmetric matrix with (i, j)-entry Qij given by (4.5), and 
where H is a (2n) x 2 matrix, in which the (i, j)-entry Hij is equal to the number of paths 
from Ai to Sj, i = 1,2, ... , 2n, j = 1, 2. (It is the negative value of the Pfaffian because of 

the sign (—1)^) on the right-hand side of (4.4), as we have q = 2.) 

In particular, using the fact that the number of lattice paths on the integer lattice Z 2 
between two given lattice points is given by a binomial coefficient, we have 

On the other hand, substituting Ai = (—i, i) and I s = (—1, s) in (4.5), we have 

l<s<t<2x+2n 

- £ ((irOCiU^C-T 

l<s<t<2x+2n V V 7 VJ 7 VJ 7 V ' 

- e (1:00:0- £ (;:0G:J 

l<s<t<2x+2n v 7 VJ 7 l<s<t<2i+2n v 7 v 

■ e C)C:0- s G)C:0 

l<t<2x+2n V 7 VJ 7 l<t<2x+2n VJ 7 V 7 

where we used the well-known identity 
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to obtain (4.9). We may obtain an alternative expression for Q i} j by replacing |(*) = j(*Z}) 

in the last expression by X Ez=o (*7 J ') G-7-i)' tms ec l uan ty being true because of the Chu- 
Vandermonde summation (cf. e.g. [16, Sec. 5.1, (5.27)]). Thus, we arrive at 

«-g'E¥f;')(.!;-,)C) 

the last line again being due to (4.11). 

To summarize, we have obtained the following result. 

Lemma 6. For all positive integers n, x and nonnegative integers k, we have 

M(F n , x \ < 2 (k)) = - Pf M n (x), (4.13) 

where M n (x) is given by (4.6), with Qij defined in (4.10) or (4.12), and Hij defined in (4.7) 
and (4.8). 

5. Proof of Theorem 4 
In the sequel, we shall interpret sums by 

n-i f Efc=i E xpr(/c) n>m 

Expr(A;) = I n = m (5.1) 

k=m { - Er=n Expr(fc) n< m. 

In particular, using this convention, the expression for Qij given in (4.10) makes sense for 
negative integers x also (in which case the upper bound in the sum can be negative) and is 
actually equal to the expression in (4.12). It is the latter fact that we shall frequently make 
use of. 

Our proof of Theorem 4 involves a sequence of five steps. By Lemma 6, we know that 
the number that we want to compute is the negative of a Pfaffian. We shall frequently use 
the fact (4.3) that the square of the Pfaffian of a skew-symmetric matrix is equal to its 
determinant. 

By its definition, Pf M n (x) is a polynomial in x. In Step 1 we prove that 

detM n (x) = detM n (-2n - x). 
With d denoting the degree of Pf M n {x) as a polynomial in x, this implies that 

Pf M n {x) = (-l) d Pf M n (-2n - x). (5.2) 
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Subsequently, in Step 2 we show that 

n 

Y\ (x + s) 2n _ 2s+1 

s=l 
s^n — k 

divides det M n (x) as a polynomial in x, while in Step 3 we show that 

n-1 

II^ + S + \)ln-2s 
s=l 

divides detM n (x). Both combined, this proves that 

n n— 1 

JJ (X + S) 2n -2s+l Y[(X + S + \)l n -2a, 

S = l 8=1 

s^n — k 

which is a polynomial of degree 

n 

J2( 4n - 4s + 1) - (2k + 1) = n(2n - 1) - (2k + 1), 

8 = 1 

divides PfM n (x) as a polynomial in x. The computation in Step 4 then shows that the 
degree of Pf M n (x), as a polynomial in x, is at most 2n 2 + n — 4k — 3. Altogether, this 
implies that 

n n 

-PfM n (x) = P„(x) JJ (x + s) 2 n-2s+lY[(x + S+±) 2 n-2s, (5.3) 

s=l s=l 
s^n— k 

where P n (x) is a polynomial in x of degree at most 

2n 2 + ?i - 4k - 3 - n(2n - 1) + (2/c + 1) = 2n - 2k - 2. 

In Step 5, we determine the value of P n (x) at x = 0, —1, . . . , — n + k + 1 (see (5.29)). 
The corresponding calculations make use of an auxiliary lemma due to Mehta and Wang 
[27], see Theorem 7 and Corollary 10 in Section 6. By (5.2), this gives us at the same 
time the value of P n (x) at x = —2n, —2n + 1, . . . , — n — k — 1. In total, these are 2n — 2k 
explicit evaluations of P n (x) at special values of x. Given the fact that the degree of P n (x) 
is at most 2n — 2k — 2, they determine P n (x) uniquely, and an explicit expression for P n (x) 
can be written down using Lagrange interpolation. If this is substituted into (5.3), then 
the evaluation of — Pf M n (x) is complete. After some manipulations, one arrives at the 
expression in (3.2). 

Step 1. detM n (x) = det M n (— 2n — x). We prove this claim by transforming, up to sign, 
M n (x) into M n (—2n — x) by a sequence of elementary row and column operations (which, of 
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course, leave the value of the determinant invariant). To be precise, for % = 2n, 2n — 1, . . . , 2 
(in this order), we add 



£(*_i)'(rowa) 

a=l N ' 



to row i, and then for j = 2n, 2n — 1, . . . , 2, we add 



Efl.J)- (column b) 

6=1 ^ ' 



to column j. Let Mn\x) denote the matrix which arises after these row and column 
operations. According to (4.9), the (i, j)-entry in Mn\x) is 



for 1 < i, j < 2n. By (4.7) and (4.8), for 1 < % < 2n and j = 2n + e, e = 1, 2, the (i, j)-entry 
of M^\x) is 

E(::i)( x ;: 2 ;-- £ 1 )' (5 - 5) 

o=l x ' x ' 

and, for 1 < j < 2n and i = 2n + e, e = 1, 2, it is 



By Chu-Vandermonde summation, we have 



, /t + 7 \ = ' /£ + 7 \ = / t + i + >y-l 

^-^ \a — I / \a + 1]J ^-^\i — a)\a J rrj) \ i + 1] 



a=l 
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whence the expression (5.4) simplifies to 




Here, we used the identity 




to obtain (5.7), and our convention (5.1) for sums to obtain (5.8). Comparison with (4.9) 
shows that this last expression is, up to the sign (— 1)* +J , exactly Qij with x replaced by 
— 2n — x. In a similar vein, the expression (5.5) simplifies to 

('T^- 2 )-^*?-*-?'*} (5 - 10) 

while expression (5.6) simplifies to the same expression with i replaced by j. Upon setting 

e = 2, this shows that the (i,2n + 2)-entry in Mn\x) is, up to the sign (—1)*, identical 
with the (i, 2n + 2)-entry in M n (—2n — x), with an analogous statement being true for the 

(2n + 2, j)-entry of M n X \x) and the (2n + 2,j)-entry of M n (-2n - x). 

We do one last row and one last column operation: in Mn (x), we add the last row to 

(2) 

the next-to-last row, and we add the last column to the next-to-last column. Let Mn (x) 

(2) 

denote the resulting matrix. By (5.10), for i = 1, 2, . . . , 2n, the (i, 2n + l)-entry of M n '(x) 
is equal to 

(-^-^ +( -K X ^-2^^-^t X ^-^- (5 ' U) 
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which is, up to the sign (— l) l+1 , exactly the (i, 2n + l)-entry in M n (— 2n — x). An analogous 

f 2) 

statement is true for the (2n + 1, j)-entries of M\ (x) and M n {— 2n — x). 

In summary, as the two by two block in the lower right corner of M n (x) consists of zeros, 
the computations (5.9)-(5.11) show that the (i,j)-entry of Mn\x) is (— 1)* +J times the 
(z, j)-entry of M n (—2n — x). This implies 

detM n (x) = detM^ 2) (x) = det M n (-2n - x), 

as claimed. 

n 

Step 2. J j (x + s) 2 n _ 2s+1 divides det M n (x) . We begin by observing that the product 

s=£n — k 

in the claim can be also rewritten as 



n n 2n—l 

II ( X + S ) 2 2n-2s+l = U ( ~ X + S ^ S+2X{S<n ~ k) ~ 2 II (x + S ) 4n - 2s+2 X (s>n+k )- 2 , (5.12) 
s=l s=l s=n+l 

s^n — k 

where x(-4) = 1 if A is true and x(^4) = otherwise. In view of Step 1, it suffices to establish 
that 



\2s+2x(s<n-k)-2 

s=l 



divides detM n (x). 

Now let s, a and 6 be integers with 1 < s < n and 1 < a < b < 2n. We claim that 



Yl ( 6 _ a )(row i of M n (-s)) = 0, (5.13) 

i=a ^ 7 

as long as 

(A) a — b < 2n — 2s < a, and 

(B) either b < 2k, or b > 2k + 3 and a-& + fc+l<ra-s<a-A;-l. 

Indeed, if we specialize (5.13) to the j-th column, where j < 2n, we obtain, using the 
expression (4.10) for Qij, 

e (!::)*, =eC::)e*^QC 

i=a V 7 i=a V 7 t=l \ / V 

2 y^ s y-v /& — a\ f ft\ ft — l\ ft — 1\ ft 



t=l i=a 
2n-2s 



b — i J \\i J \j — 1 / \i — 1 J \j 
b-a + t\ft-l\ fb-a + t-l\ ft 



few & Aj-v v 6-i / v/// (5 . 14 ) 
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Here we used Chu-Vandermonde summation (cf. [16, Sec. 5.1, (5.27)]) in the last line. Since, 
for a — b < 2n — 2s < a (which is condition (A)), the binomial coefficients containing the 
parameter b in expression (5.14) are identically zero throughout the summation range, it is 
clear that the corresponding sum vanishes. 

On the other hand, if we specialize (5.13) to the (2n + l)-st column, we obtain, again 
using Chu-Vandermonde summation, 



b 



b 



b 



i=a 



b — a\ fn — s — k — 1 



i \ i-2k-l 



b — a + n — s — k — 1 
b-2k-l 

which vanishes for b < 2k, and for b > 2k + 2 and <b — a + n — s — k — 1 < b — 2k — 1, the 
last inequality being equivalent to a — b+k+l<n— s < a — /c, and if we specialize (5.13) 
to the (2n + 2)-nd column, we obtain 



b 

E 

i=a 



b — a 
i — a 



i.2 



b 



b — a\ ( n — s — k — 1 
^ \ b-i) { i-2k-2 



— a + n — s — k — 1 
b-2k-2 

which vanishes for b < 2k — 1, and for b > 2k + 3 and < b — a + n — s — k — 1 < b — 2k — 2, 
the last inequality being equivalent to a — b + k+ l<n — s < a — k — 1. This establishes 
our claim. 

In order to prove that (x + s) 2s divides det M n (x) for 1 < s < n — k, we use (5.13) with 
a = 2n — 2s + 1 and 2n — 2s + 1 < b < 2n. It is not difficult to see that for these choices of 
s, a and b the conditions (A) and (B) are satisfied, so that we obtain 2s linear combinations 
of the rows that are linearly independent (as, for our choices of a and 6, the coefficients 
in (5.13) form a triangular array) and vanish when x = —s. This implies divisibility by 
(x + s) 2s (cf. e.g. [22, Lemma in Sec. 2]). 

To prove that (x + s) 2s ~ 2 divides detM n (x) for n — k < s < n, we use (5.13) with 
a = 2n — 2s + 1 and 2n — 2s + 1 < b < 2k on the one hand, and with a = n + k — s + 2 
and 2k + 3 < b < 2n on the other hand. Again, it is not difficult to see that for both 
types of choices of s, a and b the conditions (A) and (B) are satisfied, so that we obtain 
(2k + 2s — 2n) + (2n — 2k — 2) = 2s — 2 linear combinations of the rows that are linearly 
independent and vanish at x = — s. In the same way as before, this implies divisibility by 
(x + s) 2s ~ 2 . 

n-1 

Step 3. jfj (x + s + \)\ n _2s divides det M n (x). We begin by observing that the product 
in the claim can be also rewritten as 

n — 1 n — 1 2n — 2 

n (* + * + \)in-2s = n (* + s + i) 2s n & + * + \r- 2s - 2 - 

s=l s=l s=n 
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In view of Step 1, it suffices to establish that 

n-l 



s=l 



divides det M n {x). 

In order to prove the claim for s < n — k, we shall show that 



(n-s-k)- (row (2n - 2s + 1) of M n (-s - |)) 

1 
2 

2n-2s-l 

2 2n_2s_i + 2 



+ ^(n - s - k) ■ (row (2n - 2s) of M n (-s - §)) 

2n— 2s — 1 / x j 

+ E o2n- 2a l <+2 (row»ofM w (- 8 -l))=0, (5.15) 



i=i 
and that 

(row i of M n (- S - I)) + 21 + ^Z^Z^ ~ 3 ( row (< - !) of M ^" s " 3)) 
(2z + 2s-2n-2/c-3)(2z + 2s-2?i-2/c-5) , , , , 

+ 4( < -a-i)( < -2 t -2) 1 < row (! - 2) of M " ( - s - ^J 6 » 

for i = 2n — 2s + 2, 2n — 2s + 3, . . . , In. As these are linearly independent row combinations, 
the claim will follow. 

In order to prove the claim for s > n — k, we shall show that 

In — Is — 1 / ^ j 

E 2 2n-2s-i-i (rowzofM n (- g -i))=0, (5.17) 
i=l 

that 

(row % of M n (-s - \)) = (5.18) 
for i = 2n — 2s, 2n — 2s + 1, . . . , 2/c, and that 

(row i of M n (-s - §)) + 2Z + ^ZllZT ~ 3 ( rOW ( Z " ^ ° f M ^" S " 3» 
(2i + 2s-2n-2/c-3)(2z + 2s-2n-2A;-5) , . . „ , . 

+ ~ 4( < -a-i)( < -2 t -2) < row (I - 2) of MJ ~ S - 

for z = 2/c + 3, 2/c + 4, . . . , 2n. Again, as these are linearly independent row combinations, 
the claim will follow. 

Let first s > n — k. We start with the proof of (5.17). Specializing (5.17) to the j-th 
column, j = 1, 2, ... , 2n, by (4.9) we see that we must prove the identity 
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In order to see that this is indeed true, we first extend the sum over i to the range % 
0, 1, . . . , 2n — 2s — 1, thereby obtaining 

2n-2s-l ( — 1)* / 2n_2s_1 /j\ /-j- _ ]\ 2n— 2s — 1 _ 

E 2 2»-2.-i-l E yO-V" S Vi-00 



1 2 "^- 1 /t - 1 

o2n-2s-l Z^ \ a _ 1 



Z t=l 

2n-2s-l / lXi /2n-2s-l / , N / , N 2n-2s-l / -, \ / . 
2-^, 92n-2s-j-l I Z^ W j I i _ 1 j Z^ 

1 /2n - 2s - 1 



2 2n-2s-i-l I U/\o_iy U_l/\ ? - 



22n-2s-l I j 

for the left-hand side of (5.20). Next we interchange the sum over % with the sums over t, and 
subsequently we evaluate the (now inner) sums over % by means of the binomial theorem. 
In this manner, the left-hand side of (5.20) becomes 



1 2n-2s-l 1 2n-2s-l 

22^ E (i-^Qllj + ^izs E (i-^Q) 



1 /2n - 2s - 1 



22n-2s-l I j 

1 2 "f "'( i)' (t " 1)! (j - M) !— " 2s " ^ 

22n-2s-l Z^ V J 2 2n ~ 2s - 1 V j J 

MT'MCHi'Kr'))-* 



2 2 

as desired. 

On the other hand, specializing (5.17) to the j-th column, j = 2n + 1, 2n + 2, by (4.7) 
and (4.8) we obtain 

1« O r-, 1 . 

(-1)* ^71-)fc-S-| N 



2n— 2s — 1 / \ 4 / r S 

(-1)* fn-k- s- 6 

2 2n-2s-i-l [ j _2k-E 

i=l 



where £ = 1,2, which is indeed zero since the binomial coefficient always vanishes because 
of % < 2n — 2s — 1 < 2k + e, the last inequality being due to our assumption s > n — k. 
That (5.18) holds can be easily checked by inspection. 

For the proof of (5.19), we observe that we have Qi,j\ x= _ s _i = for all i > 2n — 2s, 

because in this case the appearance of the binomial coefficient (*) in the sum in formula 
(4.10) implies that all summands of this sum vanish. In its turn, this entails that the left- 
hand side of (5.19) specialized to the j-th column, where 1 < j < 2n, is trivially zero 
since 

% > i - 1 > % - 2 > 2k + 1 > 2n - 2s + 1 > 2n - 2s, 
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by our assumptions. To see that the left-hand side of (5.19) is as well zero when it is 
specialized to the (2n + l)-st or (2n + 2)-nd column amounts to a routine verification using 
the expressions (4.7) and (4.8) for the corresponding matrix entries. 

We now assume that s < n — k and turn our attention to (5.15). The reader should notice 
that the relations (5.15) and (5.17) are relatively similar, the essential difference being the 
two extra terms in (5.15) corresponding to the (2n — 2s)-th and the (2n — 2s + l)-st row, 
respectively. If 1 < j < 2n, the proof of relation (5.15) specialized to column j is therefore 
identical with the proof of relation (5.17) specialized to column j, because the entries in the 
first 2n columns of the (2n — 2s)-th and the (2n — 2s + l)-st row evaluated at x = —s — | are 
all zero. (The reader should recall formula (4.10).) To show the relation (5.15) specialized 
to the (2n + l)-st respectively to the (2n + 2)-nd column requires however more work. We 
have to prove 



(n 



k) 



n — s — k 
2n-2s- 2k 



3 
2 

e + 1 



If n 
+ 2\2n 



s — k 
2s - 

2n-2s-l 



3 
2 

2k -e 



+ 



E 



i-iy 



n 



k 



2 2n-2s-t+2\ i _2k- £ 



o, 



where e = 1,2, respectively, after simplification, 

(n - s - k)(e - 2) (-n + s + k + e - \)2n-2s-2k-e 



(2n-2s-2k-e + l)\ 

2n-2s-l 

+ E 



i=i 



22n-2s-i+2 



n — s — 
i - 2k 



k 



0. (5.21) 



We reverse the order of summation in the sum over i (that is, we replace i by 2n — 2s — i — 1) , 
and subsequently we write the (new) sum over i in standard hypergeometric notation 



ai, . . . , a f 



E 



b 1 ,...,b q 

Thereby we obtain 

(n - s - k)(e - 2) (-n + s + k + e - \)2n-2s-2k-e 



Z 



^ ml (6i) m • • • (b q ) m 



(2n-2s-2k-e + l)l 



(-n + k + s + e + \) 2n -2s-2k-e-i 
8(2n-2s-2fc-e-l)! : 



>.Fi 



l,-2n + 2fc + 2s + e + l. 1 
-n+k+s+e+\ '2 



(5.22) 



(5.23) 



for the left-hand side of (5.21). 

If e = 2, then the 2-F\-series in (5.23) can be evaluated using Gaufi' second 2 i 7 'i-summation 
(cf. [30, (1.7.1.9); Appendix (III.6)]) 



a,-N 



1 



2^2 2 Z 







(I) 



IFII 



N/2 



if iV is an odd nonnegative integer, 
if iV is an even nonnegative integer. 



(5.24) 



N/2 
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As a result, in this case, the expression (5.23) vanishes, whence (5.21) with e = 2 is satisfied, 
and thus relation (5.15) specialized to the (2n + 2)-nd column. 

If e = 1, the 2-Fi-series in (5.23) cannot be directly evaluated by means of GauB' formula. 
However, we may in a first stage apply the contiguous relation 



2F1 



a, b 



a, b — 1 



+ — 2^1 

c 



a + 1,6 
c+1 



to transform (5.23) into 

(n-s-k) (-n + s + h + \) 2n -2s-2k-i 



{2n~2s-2k)\ 



(-U + k + S + \)2n-2s-2k-2 

8(2n-2s-2A;-2)! 



l,-2ra + 2fc + 2s + 1 1 



-n + fc + s + | 



'2 



1 



2n-2/c-2s-3 



2^1 



2, -2n + 2k + 2s + 2 _ 1 
-n + k + s+ ^ '2 



Both 2-Fi-series in the last expression can now be evaluated by means of Gaufi' formula 
(5.24). The first series simply vanishes, while the second series evaluates to a non-zero 
expression. If this is substituted, after some simplification we obtain 



(-71 + S + k + |) 2 n-2s-2fc-l (-W + k + S + § ) 2 n-2s-2fc- 



4(2n- 2s - 2k-l)\ 



8(2n-2s -2k-2)\ 



= 0. 



This shows that for e = 1 the expression (5.23) vanishes as well, whence (5.21) with e = 1 
is satisfied, and thus also relation (5.15) specialized to the (2n + l)-st column. 

The verification of (5.16) is completely analogous to that of (5.19) and is left to the reader. 

Step 4. Pf M n (x) is a polynomial in x of degree at most 2n 2 + n — 4k — 3. By (4.12), 
Qij is a polynomial in x of degree i+j. On the other hand, the degree of H^i in x is clearly 
i — 2k — 1, while the degree of H ij2 is i — 2k — 2. It follows that, in the defining expansion 
of the determinant detM n (x), each nonzero term has degree 



2n 



2n 



2(2k + 1) - 2(2A; + 2) = An 2 + 2n - 8k - 6. 



The Pfaffian being the square root of the determinant (cf. (4.3)), the claim follows. 

Step 5. Evaluation of P n (x) at x = 0, —1, . . . , — n + k + 1. The polynomial P n (x) is 
defined by means of (5.3). So, what we would like to do is to set x = — s in (5.3), s being 
one of 0, 1, . . . , n — k — 1, evaluate Pf M n (— s), divide both sides of (5.3) by the products on 
the right-hand side of (5.3), and get the evaluation of P n (x) at x = —s. However, the first 
product on the right-hand side of (5.3) unfortunately is zero for x = —s, 1 < s < n — k — 1. 
(It is not zero for s = 0.) Therefore we have to find a way around this difficulty. 
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Fix an s with 1 < s < n — k — 1. Before setting x = — s in (5.3), we have to cancel (x + s) f 
(see (5.12)) on both sides of (5.3). That is, we should write (5.3) in the form 

Pn{x) = - ]_ Vs PfM n ( X ) 

(x + s) s 

n—k—1 n 2n—l 

x n (x+e)- £ n (x+£)- £+i n ( X +^-2»+^>n+*)+i 

1=1 £=n-k+l e=n+i 



X 



lJ(x + e+ \)£_ u , (5.25) 



and subsequently specialize x = —s. However, in order to be able to perform this step, we 
need to evaluate 

( -^—?iM n {x) 

(x + S) s 

In order to accomplish this, we apply Lemma 11 with iV = 2n + 2, a = 2n — 2s, b = 2n, and 
A = M n (x). Indeed, (x + s) is a factor of each entry in the i-th row in matrix M n (x), for 
i = 2n - 2s + 1, 2n - 2s + 2, . . . , 2n. We obtain 

1 



(x + s)' 



where 

with Q being given by 
and H by 



PfM n (x) 



Q 



-Pf(Q) Pf(5), 



Q H 
-H l 



(5.26) 



(5.27) 



Q=(Qu\ x =-.) 1 

H=(H iJ \ x _ 8 ) i 



<i,j<2n-2s 



<i<2n-2s, l<j<2 



and where 



5 = 



1 



x + s 



Q 



i+2n-2s,j+2n-2s 



l<i,j<2s 

We point out that (5.26) also holds for s = once we interpret the Pfafhan of an empty 
matrix (namely the Pfaffian of S) as 1. In particular, under that convention, the arguments 
below can be used for < s < n — k — 1, that is, including s = 0. 

We must now compute Pf(Q) and Pf(S'). We start with the evaluation of Pf(S'). It 
follows from (4.12) that the (i, -entry of S is given by 



i+2n-2s-l 

s- ■ = \ ( 3 - 1 

hJ ^ y 1 i + 2n-2s\i + 2n-2s-l-\ 

1=0 



j + 2n-2s-l \ fl+j + 2n-2s 

I 



(2n-2s + l)\(l + j 
(1 + j + 2n-2s + l)\ 
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If we write this using hypergeometric notation, we obtain the alternative expression 

' 1 - i - 2n + 2s, 1 + j + 2n - 2s, j 



s = (~ 1 ) J + 1 (J -i)i+2n-2s 

hJ (2n - 2s + j + 1)1 (j)i-j+2n-2 fl +l 



3 r 2 



l-i+j,2+j + 2n-2s 



;i 



Rewrite this expression as the limit 

s = Hm (-l) 3+1 (j -i) i+2 n-2s p ll-i~2n + 2s,l + j + 2n-2s,j _ 

<-J £ ™o (2n-2s + j + l)!(j) l - J+ 2n-2 S +i 3 1 l-z+j,2 + £ + j + 2n-2s ' 

Now we apply one of Thomae's 3 F2-transformation formulas (cf. [3, Ex. 7, p. 98]) 



3F2 

Thus, we obtain 
Sij = lim 



a, b, c 
d,e ' 



r(e)T(d + e-a-b-c) 
T(e-a)T(d + e-b-c) 



3 r 2 



a, —b + d, —c + d _ 
d, — b — c + d + e ' 



(-l) J+1 r(2n - 2s + e + 1) T(2n - 2s + j + e + 2) (j - i) l+2n _ 



2 s 



e^o T(e - i + 2) T(4n - As + i + j + e + 1) (2n - 2s + j + 1)! (j) i _ j+2n -2 S +i 

1 - i - 2n + 2s, -i - 2n + 2s, 1 - i 



1 - i + j, 2 + e - i 



x 3^2 

or, in usual sum notation, 

1 (-1) J+1 (j - i) T(2n - 2s + e + 1) T(2n - 2s + j + e + 2) 



;i 



s 



1,3 



lim > 



Z=0 



r(/ - i + e + 2) T(4n - 4s + % + j + e + 1) 

_ (1 - i)i (/ - i + j + l) 2 n-2s+i-*-i (2n - 2g + i - l)i 
/! (2n - 2s + j + 1)! (j) 2n _ 2s +j-j-/+i 

Because of the term F(l — i + e + 2) in the denominator, in the limit only the summand for 
l = i — l does not vanish. After simplification, this leads to 



{-l) i+J (j - i) (2n - 2s + i - 1)1 (2n - 2s + j - 1)1 



We must evaluate the Pfaffian 



{An -As + (2n - 2s + 1)1 



Pf (Sn). 



By factoring some terms out of rows and columns, we see that 
Pf (S l , J ) = (-iy(2n-2s + l)r s 

1<i,j<2s 



X 



2s 



H(2n-2s + i-iyA Pf 



i=i 



J ~ * 



i<i,j<2s \ (An - As + % + j)\ J ' 
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a. A lozenge tiling for the degenerate region 

Figure 7 



b. Forced lozenges in case x = 



This Pfaffian can be evaluated in closed form by Corollary 10 in the next section. The result 
is that 



Pf(S) = (-l) s (2n -2s + I)!" 5 ( JJ (2n - 2s + i - 1)! ] J 



2s 



s-1 



(2i + l) 



i=0 



(4n-2s + 2z + l)! 



. (5.28) 



We finally turn to the evaluation of det(Q). If we compare (5.27) with (4.6), then we 
see that Q = M n _ s (0). Hence, using Lemma 6 with n replaced by n — s and with x = 0, 
we see that — Pf(Q) is equal to M(F n _ S)0 \ <2(^))- (The reader should recall the definitions 
of the region F njX and of the triangular hole <2(k) given in the introduction, see again 
Figure 4.) Figure 7. a shows a typical example where n — s = 5 and k = 2. Since this region 
is degenerate, there are many forced lozenges, see Figure 7.b. The enumeration problem 
therefore reduces to the problem of determining the number of symmetric lozenge tilings of 
a hexagon with side lengths 2k, 2k, 2, 2k, 2k, 2. This number is given by formula (3.1) with 
n = k and x = 1. If we substitute this in (5.26), together with the evaluation (5.28), then, 
after some manipulation, we obtain 



(x + s)< 



PiM n (x) 



(-1) 



4k + 1 



(2s) 



2k J (2n - 2s + l)\ s 2 s s\ 

2s v , s — 1 

x ( n(2"-2*+i-i)!)(n 



(2i)! 



LA (4n-2s + 2i + l)! 

i=0 
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Hence, by inserting this in (5.25), we have 
'4fc + l\ (2s)! 



P n (-a) = (-iy 



x 



2k J (2n- 2s + l)! s 2 s s! 

2 * s x ' {2i)\ 



Yl(2n-2s + i-l)\ J] 



. . n (4n-2s + 2i + l)! 



n—k—l n 2n—l 

X 



J] (s + £)- £ Jl (-S + £)- £+1 J]' (_ s + £)-2n+^>n+fc) + l 



1=1 e=n-k+i e=n+i 



n 

X 



l[(- s + £ + (5.29) 



This completes the proof of Theorem 4. 



6. AN AUXILIARY DETERMINANT EVALUATION, 
AND AN AUXILIARY PFAFFIAN FACTORIZATION 

Mehta and Wang proved the following determinant evaluation in [27]. (There is a typo 
in the formula stated in [27, Eq. (7)] in that the binomial coefficient is missing there.) 

Theorem 7. For all real numbers a, b and positive integers n, we have 
det ((a + j-z)r(6 + z+j)) 

0<i,j<n — 1 

(n-1 \ n / \ 

Jl i! r(6 + i) J £(-l) fc Q ((6 - a)/2) fc ((6 + a)/2) n _ fc , (6.1) 
i — Iz := 

as /ong as i/ie arguments occurring in the gamma functions avoid their singularities. 

The sum on the right-hand side of (6.1) can be alternatively expressed as the coefficient 
of z n /n\ in 

(l + z )(«-*)/2(l_ z )(-«-«0/2. 

Therefore, in the case a = we obtain the following simpler determinant evaluation. 
Corollary 8. For all real numbers b and positive integers n, we have 

n-l 



det ((j - i) T(b + i + j)) = x{n is even) ( [ J i\ T{b + i) J 



(n/2)! 



as /ono as £ae arguments occurring in the gamma functions avoid their singularities. Here, 
as before, x(^4) = 1 if A is true and x(-A) = otherwise. 

One can obtain the following slightly (but, for our purposes, essentially) stronger state- 
ment. It is stated as Eq. (4) in [27], with the argument how to obtain it hinted at at the 
bottom of page 231 of [27]. Since, from there, it is not completely obvious how to actually 
complete the argument, we provide a proof. 
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Proposition 9. For all real numbers b and positive even integers n, we have 

— — 1 

Pf (U-i)T(b + i + j)) = IJ(2i + l)!r(6 + 2i + l), (6.2) 

0<i,j<n— 1 

i=0 

as long as the arguments occurring in the gamma functions avoid their singularities. 

Proof. Since the Pfaffian of a skew-symmetric matrix equals the square root of its determi- 
nant (cf. (4.3)), the formula given by Theorem 8 yields, after a little manipulation, that 

Pf Uj -i)T(b + i+j)) =e FT (2z + l)!r(6 + 2z + l), (6.3) 

,7<n— 1 - LJ - 



0<i,j<n- 

^=0 

where e = +1 or e = —1. In order to determine the sign e, we argue by induction on (even) 
n. Let us suppose that we have already proved (6.2) up to n — 2. We now multiply both 
sides of (6.3) by b+ 1 and then let b tend to — 1. Thus, on the right-hand side we obtain the 
expression 

^ i=0 ' ^ z=l ' 

On the other hand, by the definition of the Pfaffian, on the left-hand side we obtain 

E sgn Tallin | (6+1) J] (j-i) T (b + i + j)\ (6.5) 

TrEM[0,...,n-l] V i<j J 

i,j matched in tt 

(with the obvious meaning of M. [0, . . . , n — 1] ; cf. the sentence containing (4.2)). In this sum, 
matchings n for which all matched pairs i,j satisfy i +j > 1 do not contribute, because the 
corresponding summands vanish. However, there is only one possible pair i,j with < i < j 
for which i + j < 1, namely (i, j) = (0, 1). Therefore, the sum in (6.5) reduces to 

E sgn7r , ( 6 lim i (6+l)(l-0)r(6+l)) (j ~ i)T(i + j - I) 

■K'€M[2,...,n-l] ' i<j 

i,j matched in ir' 

= Pf Oj-i)T(i + j-l)) 

2<i,j<n — 1 

= n< Pf (U-i)T(i+j + S)), 

0<i,j<n — 3 

where the next-to-last equality holds by the definition (4.2) of the Pfaffian. Now we can use 
the induction hypothesis to evaluate the last Pfaffian. Comparison with (6.4) yields that 

£ = +1. □ 

By using the reflection formula (cf. [2, Theorem 1.2.1]) 

T(x)T(l-x) = 

Sin 7TX 

for the gamma function, and the substitutions % — > n — i — 1 and j — > n — j — 1, it is not 
difficult to see that Proposition 9 is equivalent to the following. 
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Corollary 10. For all positive even integers n, we have 



Pf 



J - * 



o<i,j<n-i \T(b + i + 



n 

i=0 



(2i + l)! 



T(b + n + 2i- 1) 



We close this section by proving a factorization of a certain specialization of a Pfaffian 
that we need in Step 5 in Section 5. 

Lemma 11. Let N, a, b be positive integers with a < b < N , where N and b — a are even. 
Let A = (Aij)i<ij<N be a skew- symmetric matrix with the following properties: 

(1) The entries of A are polynomials in x. 

(2) The entries in rows a + 1, a + 2, . . . , b (and, hence, also in the corresponding columns) 
are divisible by x + s. 

Then 

1 



(x + s)( b " a )/ 2 



Pf A 



Pf A-PiS, 



(6.6) 



where A is the matrix which arises from A by deleting rows and columns a + 1, a + 2,. ..,6 
and subsequently specializing x = —s, and 



5 = 



1 



x + s 



a+l<i,j <b 



Proof. By the definition (4.2) of the Pfaffian, we have 



(x + sY h - a )/ 2 



Pf A 



J2 sgnyr J] 



( X + 8 )(b-a)/2 ^ 11 ~*>3 

n€M[l,...,N] i<j 

i,j matched in tt 



Let A4i denote the subset of M. [1, . . . , N] consisting of those matchings that pair all the 
elements from {a + 1, a + 2, . . . , b} among themselves (and, hence, all the elements of the 
complement {1, 2, . . . , a, 6 + 1, 6 + 2, . . . , iV} among themselves). Let M.2 be the complement 
M[l,...,N]\Mi. Then 



(X + S)( & " a )/ 2 



Pf A 



( x + 8 )V«)/2 £ sgn7r n — 

. ' nEMi i<j 

i,j matched in tt 

+ f {x + 8) Va)/2 E ^ n a ^ 

\ v 7 7i-eX 2 i<j 

i,j matched in 7r 



(6.7) 



Each term in the sum in the second line of (6.7) vanishes, since the product contains more 
than (6 — a)/2 factors Aij that are divisible by x + s. On the other hand, every matching 
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7T in Mi is the disjoint union of a matching it' e M[l, 2, . . . , a, 6 + 1, 6 + 2, . . . , N] and a 
matching f" G Af[o+l,o + 2,...,6]. If we also use the simple fact that sgn7r = sgn n' ■ sgn n" 
(as there are no crossings between paired elements of it' and paired elements of 7r"), then 
we obtain 



(x + s)( b " a )/ 2 



Pf A 



i y 

(x + s)( b - a y 2 ^ 

V ' n'eM[l,...,a,b+l 

n"£M[a+l,...,b] 



n ^ 



sgn tt' ■ sgn 7r" 



n a > 



1,3 I \ 11 ' J 

i,j matched in tt' i,j matched in tt 



sgnyr' Y[ Aij 

7r'eM[l,...,a,b+l,...,N] i<j 

i,j matched in tt' 



e n 

7r"€M[o+l,...,6] i<j 

i,j matched in it' 



X + S 



By the definition (4.2) of the Pfafnan, the last expression is exactly the right-hand side of 
(6.6). □ 

7. Proofs of Theorems 1 and 2 
In our proofs we make use of the following lemmas. 

Lemma 12. Let (3 be a real number with either (3 > or j3 < — 1. Then, for all sequences 
(Pn)n>i with j3 n — > j3 as n — > oo, we have 



lim -j= 5-F4 



— 2n, 2 , — 7i + + 1, — n + + 1, /3 n n 
— 2n + i, — n — fc, — n — k, (3 n n + 1 



1 V2(l-a) 4k + 2 
+ v/7ra(2-a) 



1 = / ~ . - , ^ 

(7.1) 

where, on the left-hand side, we used again the standard notation (5.22) for hypergeometric 
series. 

Proof. We write the 5-F4-series in (7.1) explicitly as a sum over /: 



-k-l 



T(2n + 1) T(Z + \) T(2n -l + \)T{n- kf T(n + k-l + if 



^ r(2n - / + 1) T(i) T(/ + l) r(2n + \) T(n - k - I) 2 T(n + k + l) 2 (p n n + I) ' 



(7.7) 
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Let us denote the summand in this sum by F(n, I). We have 

^F(n, I) = F(n, I) (^(l + \) - ^{l + 1) + ^{2n - I + 1) - ^(2n - I + \) 

+ 2^(n -k-l) - 2^{n + k — I + 1) — 1 ) , 

where ip(x) := (-^T(x))/T(x) is the digamma function. Since ip(x) is a monotone increasing, 
concave function for x > (this follows e.g. from [2, Eq. (1.2.14)]), we have 

V>(/ + 1) - */>(l +^)> ^(2n - / + 1) - ifj(2n - I + \) 

for < / < n. Moreover, because of the equality ip(x + 1) = ip(x) + - (cf. [2, Eq. (1.2.15) 
with n = 1]), and since either f3 > or f3 < — 1, for large enough n we have 

ip(n + k - I + 1) > ip(n - k - I) H > ip(n - k - I) - ' 



n + k — I (3 n n + I 

Altogether, this implies that JyF(n, Z) < for < I < n — k — 1, that is, for fixed large 
enough n, the summand -F(n, /) is monotone decreasing as a function in /. In particular, for 
< I < n — k — lwe have 

< F(n,l) < F(n,0) = 1. (7.8) 
The sum (7.7) may therefore be approximated by an integral: 

n-k — l |_ lo g n J — 1 n — k— |_lognJ -1 n -k — l 

F( n ,i)= Yl F ^ l ^+ E F M + E f m 

Z=0 Z=0 Z=|k>gnJ l=n-k-\_\ogn\ 

r-n — k — [log nj —1 

= 0(logn)+ / F(n,l)dl. (7.9) 

«/ |_i°g "J —i 

The next step is to apply Stirling's approximation 

log T(z) = - log(z) - z + \ log(27r) + O (7-10) 
for the gamma function, in the form 

log T(an + bl + c)= (an + bl + c-^j (log (a + b±) + log(n) + log (l + ) 

- (an + bl + c) + J log(2yr) + O 



2 \ an + bl 

an + bl + c-^j (log(a + &£) + log(n)) 

- (an + 6/) + J log(2yr) + O 1 



2 V on + bl 
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where a, b, c are real numbers with a > 0. If this is used in the defining expression for F(n, I), 
then after cancellations we obtain 

logF(n, /) = Jlog(2) + (4k + 2) log f 1 - - ) - ^log f 2 - - 



nil \ n 



7: log ( - J - Jlog(n) - Jlog(yr) - log ( 1 + 



log 




2 \nj 2 ov ' 2 x ' \ (3 n n 

1 



+ 0\jJ+0[- _)+0[__ r )+o 

V2(l-M 4k + 2 



n 



.^(1 + ^)^(2-^), 
as long as logn < / < n — k — logn. Substitution of this approximation in (7.9) yields 



n—k — l 




(3 n n 

or, after the substitution / = an, 

n—k — l 




1 + +0(logn), 



z=o 



(p(n—K— |_ic 
/ 



(n-fc-|k>gnj-l)/n (1 - a ) 4fc + 2 



)/n (1 + fa) ^/7ra(2-a) / V ' V>g™ 




l + O +0(logn 



The assertion of the lemma follows now immediately. □ 

With no extra work we can now get an exact formula for a generalization cu/(/c;£) of 
the correlation ujf(k) described in Section 2. For any real number £ > 0, define w/(/c;£) in 
analogy to (2.1) by 

^(fcO^iimM^WM, (7 . n) 

where (£ n ) n >i is a suitable sequence of rational numbers approaching £. ("Suitable" here 
means that we have to choose £ n in such a way that £ n n is integral.) The number tvf(k; £) is 
the correlation of the triangular gap <2(k) with the free boundary, obtained when the large 
regions used in the definition are the left halfs of hexagons that are not necessarily regular, 
but have their left vertical side £ times as long as the two oblique sides. Note that, by the 
results of [9], we should expect distorted dimer statistics around the gap for £ 7^ 1. As 
Theorem 15 shows, the distortion is quite radical: 0Jf(k; £) turns out to decay exponentially 
to or blow up exponentially, according as £ > 1 or £ < 1; see also Remark 1. 
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Lemma 13. For any £ > and < k E Z, we have 
1 fAk + l\ 1 



1 A 



2k J (l + £)4fc+2 v /2+^ 



, ^ f 1 (l-«) 4fc+2 , , f 1 (l-«) 4fc + 2 
x (£ + 2) / ; da-£ i ' da 

J Jo (1 + f) V«(2-«) 7o (l-^)v / ^(2^y 

_ 1 /4fc + 1\ 1 r 1 2(1- g) 4fc + 3 

~n{ 2k J (l + 4k+2 V2 + ZJo (l+f)(l-^)v^(2^) ( ' J 

Proof. By Theorem 4 and formula (3.1), the ratio between M(F n>x \<2(&)) an d M(F n)X ) is, 
when written in hypergeometric notation, 

'4k + l\ {n + k)\{\) 2n {x + l) n _ fc _i (x + n + k + l) n -fc-i 



2/c y (x + n - fe) 2 fc+i (a; + §) 2 „ (n-k- l)! 2 (n + fc + l) n -k 

-2n, |, — n + k + 1, — n + + 1, a; 



x (x + 2n) 5F4 



271 + 1, —?i — k, —n — k, x + 1 ' 



-2n, i, — n + + 1, — n + k + 1, — 2n — a; 
— 2n + |, —?7. — fc, — n — k, —2n — x + 1 ' 



We now substitute x = £ n n in this expression. Use of Lemma 12 (which applies, as £ > 0), 
together with Stirling's formula (7.10), yields the assertion. □ 

Lemma 14. For any (3 7^ we have 

(i- a ) 4fc + 2 =da A ,^ oo _ (7J3) 

(1 + § ) V«(2 - a) VSfc 

Proof. Let J/3 (A;) be the integral on the left hand side of (7.13). The asymptotics of Ip(k) 
as k — > 00 can be readily found using Laplace's method as presented for instance in [28]. 
Conditions (i)-(v) of [28, pp. 121-122] are readily checked. By [28, Theorem 6.1, p. 125], 

the large z asymptotics of e~ zp ^q(t) dt is determined by the quantities A, /j, po and qo 
in the series expansions 

p(t) - p(a) =p (t- aY +pi(t- ay +1 + ■■■ 

and 

q(t) = q (t - a) x + qi(t - a) x+1 + ■■■ . 
Namely, under the above assumptions one has 

e-<«> f\-^ q(t)dt = T (h)^g!l + o( * V (7.14) 



In the case of Ip(k) we have p(t) = — ln(l — £), q(t) = 1 . , a = 0, and 6=1. 

( 1 - i //3)v*( 2_t ) 

These yield parameters A = 1/2, /j, = 1, p = 1, and <j» = l/v 7 ^- In addition, p(a) = 0. As 
in our case 2; = 4/c + 2, under these specializations (7.14) becomes (7.13). □ 
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Theorem 15. As k — > oo, the correlation u)f(k;£) is asymptotically 



1 1 - ^ ;A 



7r(l+0 2 v^(2TO * Vl + € 

Proof. Combine Lemmas 13 and 14 with Stirling's approximation for the binomial coefficient 

(it 1 ) in (7.12). □ 

Proof of Theorem 1. Set £ = 1 in Theorem 15. □ 
Proof of Theorem 2. Set 

£> fc :=3Ji(A0-J_ 3 (A0, 

where Ip(k) denotes the integral on the left hand side of (7.13). Recalling that ujf(k) is the 
£ = 1 specialization of Uf(k;£), we have by Lemma 13 that 

(4£; + 3)(4£; + 5) 



Uf {k+l)-u; f {k) = -^j={ 2k ) 



4(2k + 2)(2k + 3) 



1 



D k+1 + (D fc+ i - D k ) 



and thus 



^/(fc + 1) -ujf(k) 
ujf(k) 



(4fc + 3)(4fc + 5) _ ' 
4(2fc + 2)(2Jfe + 3) 



By two applications of Lemma 14 it follows that 



"dT + a ' (7 ' 16) 



*-oo. (7.16) 

Thus Dk+i/Dk — > 1 as fe — > oo, and elementary arithmetics implies that the first term on 
the right hand side of (7.15) is asymptotically — l/(2fc) as — > oo. 

To determine the asymptotics of the second term, write by Lemma 13 

D k+1 -D k = 3 [h(k + 1) - h(k)} - [I_ 3 (k + 1) - 7_ 3 (fc)] . (7.17) 

As Ip{k) is the integral on the left hand side of (7.13), we have 

r 1 (1 - a) 4k+2 

Ip{k + 1) - Ip(k) = / 1 J — [(1 - a) 4 - 1] da. (7.18) 

■/<> (i + fjv^^y 

The asymptotics of the integral in (7.18) follows by Laplace's method, in the same manner 
as the proof of Lemma 14. In this case A = 3/2, /j, = 1, and equations (7.14) and (7.18) 
impliy that 



Ip(k+l)-Ip{k)~-jgL-, k^oo. (7.19) 

Equations (7.17) and (7.19) determine the asymptotics of D k +i — D^, and combining this 
with the asymptotics of D k given by (7.16) we obtain that the second term on the right hand 
side of (7.15) has asymptotics — l/(2k) as k — > oo. The two terms on the right hand side of 
(7.15) thus have a sum that is asymptotically — l/(2fc) — l/(2fc) = —1/k, and Theorem 2 is 
proved. □ 
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